On global Hessenberg based methods for solving Sylvester matrix equations

In the first part of this paper, we investigate the use of Hessenberg-based methods for solving the Sylvester matrix equation $AX+XB=C$. To achieve this goal, the Sylvester form of the global generalized Hessenberg process is presented. Using this process, different methods based on a Petrov–Galerkin or on a minimal norm condition are derived. In the second part, we focus on the SGl-CMRH method which is based on the Sylvester form of the Hessenberg process with pivoting strategy combined with a minimal norm condition. In order to accelerate the SGl-CMRH method, a preconditioned framework of this method is also considered. It includes both fixed and flexible variants of the SGl-CMRH method. Moreover, the connection between the flexible preconditioned SGl-CMRH method and the fixed one is studied and some upper bounds for the residual norm are obtained. In particular, application of the obtained theoretical results is investigated for the special case of solving linear systems of equations with several right-hand sides. Finally, some numerical experiments are given in order to evaluate the effectiveness of the proposed methods.     

A Reliable Stochastic Numerical Analysis for Typhoid Fever Incorporating With Protection Against Infection

In this paper, a reliable stochastic numerical analysis for typhoid fever incorporating with protection against infection has been considered. We have compared the solutions of stochastic and deterministic typhoid fever model. It has been shown that the stochastic typhoid fever model is more realistic as compared to the deterministic typhoid fever model. The effect of threshold number T* hold in stochastic typhoid fever model. The proposed framework of the stochastic non-standard finite difference scheme (SNSFD) preserves all dynamical properties like positivity, bounded-ness and dynamical consistency defined by Mickens, R. E. The stochastic numerical simulation of the model showed that increase in protection leads to low disease prevalence in a population.     

Multiphase topology optimization with a single variable using the phase-field design method

In this study, a multimaterial topology optimization method using a single variable is proposed by combining the solid isotropic material with penalization method and the reaction-diffusion equation. Unlike ordinary multimaterial optimization, which requires several variables depending on the number of material types, this method intends to represent various materials as one variable. The proposed method combines two special functions in the sensitivity analysis of the objective function to converge the design variable into prespecified density values defined for each of the multimaterials. The composition constraint based on a normal distribution function is also introduced to estimate the distribution of each target density value in a single variable. It enables density exchange between multiple materials by increasing or decreasing the amount of a specific material. The proposed method is applied to structural and electromagnetic problems to verify its effectiveness, and its usefulness is also confirmed from the viewpoint of cost and computation time.     

A CCD-ADI method for two-dimensional linear and nonlinear hyperbolic telegraph equations with variable coefficients

In this paper, a combined compact finite difference method (CCD) together with alternating direction implicit (ADI) scheme is developed for two-dimensional linear and nonlinear hyperbolic telegraph equations with variable coefficients. The proposed CCD-ADI method is second-order accurate in time variable and sixth-order accurate in space variable. For the linear hyperbolic equation, the CCD-ADI method is shown to be unconditionally stable by using the Von Neumann stability analysis. Numerical results for both linear and nonlinear hyperbolic equations are presented to illustrate the high accuracy of the proposed method.     

Global weak solutions to a spatio-temporal fractional Landau–Lifshitz–Bloch equation

The present paper deals with global existence of weak solutions of a time-space fractional Landau–Lifshitz–Bloch equation involving the weak Caputo derivative and a fractional Laplacian. We use Faedo–Galerkin method with some commutator estimates in order to prove global existence of weak solutions for the model. The uniqueness is also discussed in a special one dimensional case.     

Deposition of UF5 on the surface of steel material having oxide coating in contact with gaseous UF6

When UF6 is handled in leak tight system, the amount of uranium compound formed on the surface of containers, valves, and others in UF6 handling facilities is not so significant for a short term that special attention has not been given to this problem. The present work was done to throw some more light on this problem based on the recent experiment. We discuss the possibility that the intermolecular transfer of a fluorine atom from UF6 to UF5 may participate in the formation of uranium compound. The discussion includes also the unique features contained in the experimental result, the reaction processes assumed in this problem, and the derivation of a rate equation for expressing the deposition of uranium compound. Furthermore, we propose a new method for determining nonlinear parameters included in a governing differential equation having two variables for expressing the deposition of uranium compound from experimental raw data.     

On the numerical solution of Fredholm integral equations utilizing the local radial basis function method

The current investigation describes a computational technique to solve one- and two-dimensional Fredholm integral equations of the second kind. The method estimates the solution using the discrete collocation method by combining locally supported radial basis functions (RBFs) constructed on a small set of nodes instead of all points over the analysed domain. In this work, we employ the Gauss–Legendre integration rule on the influence domains of shape functions to approximate the local integrals appearing in the method. In comparison with the globally supported RBFs for solving integral equations, the proposed method is stable and uses much less computer memory. The scheme does not require any cell structures, so it is meshless. We also obtain the error analysis of the proposed method and demonstrate that the convergence rate of the approach is high. Illustrative examples clearly show the reliability and efficiency of the new method.     

一类Hilfer型分数阶微分方程解的存在和唯一性

用Banach不动点定理和上下解方法研究了α(1<α <2)阶数和β(0≤β ≤1)类型Hilfer分数阶微分方程的解,给出了方程解的存在和唯一性,并通过例证验证了本文所得结果的有效性.